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factor needed to match maximum thickness for pointed tail bodies using 
modified slender body theory 
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body length 
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radial coordinate, (z 2 + y 2 ) 1 / 2 
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t maximum thickness of body (measured from body centerline) 

U q free- stream velocity 

U r radial velocity 

(U ,U ,U ) velocity components in (x,y,z) coordinates 
x y z 

XE amount of extension of body length 

(XpjY^jZp) coordinate system with origin at c.g, 

(XgjYgjZg) coordinate system with origin at hub position 
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interference velocity 
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lateral cyclic control angle 
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SUMMARY 


A computationally efficient body analysis designed to couple with a comprehen- 
sive helicopter analysis is developed in order to calculate the body-induced aerody- 
namic effects on rotor performance and loads. A modified slender body theory is used 
as the body model. With the objective of demonstrating the accuracy, efficiency, and 
application of the method, the analysis at this stage is restricted to axisymmetric 
bodies at zero angle of attack. By comparing with results from an exact analysis for 
simple body shapes, it is found that the modified slender body theory provides an 
accurate potential flow solution for moderately thick bodies, with only a 10%-20% 
increase in computational effort over that of an isolated rotor analysis. The compu- 
tational ease of this method provides a means for routine assessment of body-induced 
effects on a rotor. Results are given for several configurations that typify those 
being used in the Ames 40- by 80-Foot Wind Tunnel and in the rotor-body aerodynamic 
interference tests being conducted at Ames. A rotor-hybrid airship configuration is 
also analyzed. 


INTRODUCTION 


The aerodynamic interaction between the rotor and the fuselage of a helicopter 
requires careful consideration when assessing performance and loads. This interac- 
tion has been shown experimentally to have significant effects (e.g., in refs. 1, 
and 2). Ideally, analytical models of the rotor-induced effects at the body and the 
body-induced effects at the rotor need to be combined. Because of the complex nature 
of the flow field of a helicopter, however, the end result of such a combination 
would require the coupling of complex body and rotor computer codes in order to cal- 
culate the system behavior. Typical past efforts (refs. 3 and 4) have coupled poten- 
tial flow panel methods for fuselage aerodynamic analysis with simplified rotor 
analytical models. For investigations emphasizing the rotor, an efficient analytical 
model for the fuselage is desirable, rather than a complex and time-consuming paneling 
code. 


This paper presents the initial results of an investigation intended to develop 
an efficient body model for use in comprehensive helicopter analyses, such as CAMRAD 
(ref. 5), in order to calculate body-induced effects on rotor performance and loads. 
The. requirement for efficiency does not, however, diminish to any degree the need to 
maintain the accuracy of the body flow field. Computational efficiency is usually 
achieved by using simple methods; in this case, a direct rather than an inverse or 
iterative method is desired. Slender body theory is the logical starting point for 
the analysis, but modifications to the theory will be required to achieve the accu- 
racy desired for moderately thick bodies. 

This investigation will ultimately lead to a general model for the body shape, 
and will include flows at nonzero angle of attack. The' results of the present paper 
are the initial steps of the investigation, including (1) development of a model for 
a body of revolution at zero angle of attack; (2) demonstration of the accuracy and 



limits of the model; and (3) calculations of the isolated body-induced effects on the 
rotor. 

Ames Research Center is engaged in an experimental program to measure the aero- 
dynamic interaction between the rotor and fuselage. Several small-scale tests have 
been completed (refs. 6-8), and more are planned; full-scale tests will be conducted 
in the near future. For the baseline fuselage bodies, these investigations use axi- 
symmetric shapes. Hence, although the present analytical investigation is restricted 
to axisymmetric shapes for the purpose of demonstrating its accuracy and limits, the 
method is directly applicable to the configurations of this experimental program. 

The two test modules routinely used for full-scale rotor tests in the Ames 40- by 
80-Foot Wind Tunnel — the Rotor Test Apparatus (RTA) and the Easter Egg (EE) — pro- 
vide body shapes for calculation of rotor performance and loads. Calculations can 
also be made with this body model to examine the ability to obtain isolated rotor 
behavior from measurements on these standard test modules. 

This report presents the development of a modified slender body theory to calcu- 
late the flow field of an axisymmetric body at zero angle of attack. This analysis 
is coupled with a comprehensive helicopter analysis for calculation of the body- 
induced effects on the rotor performance and loads. The approach is verified by com- 
paring the results for performance and loads from the modified slender body theory 
with results from an exact analysis for simple shapes. This comparison also provides 
a first look at the source of the rotor performance and load changesi ' Next, the 
body-induced effects on the rotor are calculated for several realistic cases: a 

typical full-scale rotor test configuration, configurations representing the model and 
full-scale aerodynamic interference tests at Ames, and a rotor-hybrid airship config- 
uration. Finally, the efficiency of the modified slender^ body theory computations 
is discussed. 


BACKGROUND OF ROTOR-BODY INTERACTION 


Body-induced effects on the rotor have received somewhat less attention than 
rotor-induced effects on the body, although the entire topic of rotor-body aerody- 
namic interaction has been subject to intensive investigation only since the lat^e 
1970s. The following is a brief summary of some of the work dealing with body- 
induced effects on rotor behavior. 

Reference 2 analyzes the effects of a model Lynx fuselage on rotor behavior. 
Calculations were made of the fuselage upwash effects on blade lift, torque, and 
bending moments. The change in blade angle of attack caused by the body, as well as 
the effects of rotor-body separation on upwash velocity, is presented. Body-induced 
effects on hub forces and moments are also shown. Reference 3 shows measurements of 
the effects of fuselage width and rotor-body separation on blade moments. A model 
fuselage and a fully articulated model rotor were used. Reference 4 presents calcu- 
lations showing the effect of the fuselage on rotor blade angle of attack for two 
advance ratios. A sparse contour plot of the fuselage upwash velocity in the rotor 
plane is also shown. In addition, reference 4 provides a summary of some of the 
analytical work done in the rotor body interaction area. 

Reference 6 presents experimentally the effects of a body of revolution on rotor 
performance. Advance ratio, tip-path-plane angle, body angle, rotor-body separation, 
and hub position are varied to cover a range of configurations. In reference 7, 
however, it was concluded that the data obtained in reference 6 were not accurate 
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enough for making conclusions about the influence of the body in rotor performance. 
Reference 8 is part of the interactional aerodynamics program being conducted at 
Ames, which includes references 6 and 7. Of particular interest to the present paper 
are the results shown in reference 8 for the surface pressure distribution of an 
isolated body of revolution. 

Reference 9 shows the calculated upwash velocity contour in the rotor plane for 
a general fuselage model for several rotor-body separation distances. Also shown is 
the calculated influence of the fuselage on the local blade angle of attack and local 
blade lift for several rotor-body separation distances. Reference 10 shows the cal^ 
culations for the fuselage upwash velocities in the rotor plane and the change in the 
blade angle of attack. Two analytical models were used in calculating the change in 
angle of attack owing to the body. The more refined model showed changes of smaller 
magnitude than the less complex model. 

Reference 11 provides calculations of the effect of the fuselage on the edgewise 
bending moments for several flight conditions. The change in the local blade angle 
of attack caused by two body shapes (a helicopter fuselage and a wind tunnel test 
module) is shown. Reference 12 presents data from a model YUH-61A, The effects on 
blade bending moments and blade torsion resulting from variations in rotor-body sepa- 
ration and forward speed are presented. From the tabulated data presented in refer- 
ence 13, the effect of two model fuselages on the rotor drag-to-thrust ratio can be 
found (effects shown in ref. 10). . 

Overall, there appears to be almost no performance measurements and no data for 
axisymmetric bodies when considering body-induced effects on the rotor. There has 
been little analysis of axisymmetric shapes (ref. 11) and little use of slender body 
theory with rotor calculations. 


DEVELOPMENT OF THEORY 


In developing an analytical model for the body flow field, the primary consider- 
ation was the computational efficiency of the method. The calculation of isolated 
rotor performance and loads is already a computationally expensive task. If the com- 
bined rotor-body analysis is to be predicted for routine use, it must require the 
same order of computer time as the isolated rotor problem. A number of panel methods 
have been developed for calculating the flow field of arbitrary bodies. These methods 
typically require several times the computer time of the rotor analysis; hence their 
usefulness is limited. A more efficient method is thus required. Considering the 
sources of the large computation time required for the panel codes, it is concluded 
that an efficient technique must use on the order of 100 singularities, and must 
obtain the singularity strength directly from the body shape (without inversion or 
iteration steps). Slender body theory, using small integration steps to evaluate the 
body-induced velocities, satisfies these requirements. For the cases considered here 
(axisymmetric bodies at zero angle of attack) it is possible to develop a surface 
singularity method that would be reasonably efficient. Unlike slender body theory, 
however, such a technique would not retain its efficiency when extended to the more 
general prob lem . 
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Slender Body Theory 


The basis for the body model is slender body theory; potential flow is assumed. 
The body thickness is assumed to be much less than the body length. Since axisym- 
metric bodies are assumed, the potential and stream function of an axial source 
distribution are used (ref. 14): 


<j> (x,r) 


x f q(Odg 

4lV J A [(x - ?) 2 + r 2 ] 1/2 


(x , r ) 



q(Q (x- ?)d? 

[(* - S ) 2 + r 2 ] 1/2 


The coordinate system and geometry used are shown in figure 1. Here, q represents 
the source distribution and A and B represent the endpoints of the body. The prob- 
lem is to find a distribution q that will model the body. The following conditions 
must also be satisfied: 


= 0 (velocity normal to the body surface is zero) 
an 


(i) 



q (x)dx 


= 0 


(net source strength is zero for a closed body) 


( 2 ) 


Applying conditions (1) and (2), in addition to asstiming r/x « 1, the following 
result is obtained (see ref. 14) for details of derivation): 


q(0 = U o A' (?) 


(3) 


Here, A is the cross-sectional area of the body. For a body of revolution, 

A' (?) = 2iTr 0 (?)rQ(?) . Substituting for q in the expressions for <{> and \jj yields: 


<j> (x, r) 


'P (x,r) 



r o (?) r c>(?)d? 
[(x - ?) 2 + r 2 ] 



r o (?)r;(?)(x - ?)d? 

t(x - ?) 2 + r 2 ] 1/2 


(5) 


Note that the above two equations now include the free-stream velocity U 0 . 

To obtain the axial and radial velocity components induced by the body, the 
Cauchy- Reimann relations are used, that is. 


ai.JLii.u 

3x 2irr 9r x 
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M _ 1 M = n 

9r 2*jrr 9x r 


Since the zero streamline, ^ = 0, defines the body surface, the stream function 
rather than the potential is used in subsequent equations and calculations. Applying 
the Cauchy-Riemann relations yields 


U 

x 



r 0 (S)r’a)(x - OdK 

t(x - O 2 + r 2 ] 3 /Z 


u 


r 



r o (Or^te)r d.g 

[(X - S) 2 + r 2 ] 3/2 


Ideally, the stagnation points of the flow field are located at the nose and 
tail of the body. For bodies with blunt noses and tails, the stagnation points will 
be forward of the nose and aft of the tail, if the source distribution is allowed to 
run from the nose to the tail. The stagnation point is found by setting the equation 
for U x equal to zero at r = 0 and solving the resulting equation for x by 
iteration. Figure 2(a) shows the amount of overshoot of the stagnation point versus 
the maximum thickness of several ellipses. To correct for this o vet shoot , the limits 
of integration were modified. For bodies with blunt ends, reference 15 shows that 
good approximations for the locations of the limits of integration are points that 
are halfway between the nose and the center of curvature of the nose, and halfway 
between the tail and the center of curvature of the tail. 


Location of Integration Limits 


The required offset for the limits of integration can be derived in the follow- 
ing way. Assume some offset e at the nose. Physically, U x = 0 at the actual body 
nose, that is, U x = 0 at (x,r) = (A - e,0) . So, 

1 r B 

0 = 1 -J J r o (Or^a)a - A + e)d? 

•'A 

First consider the case of a blunt or rounded nose. For this case, r Q (A) = 0 
and r^CA) = °°, but r 0 (A)ro(a) is finite. In fact, r 0 (A)ro(A) = RCN, where vRCN is 
the radius of curvature of the nose. This can be shown by examining the equation for 
a circle with radius of curvature RCN and tangent to the body nose: 

RCN 2 = (x - RCN) 2 + r 2 (x) 


The derivative is 


0 = 2 (x - RCN) + 2r Q (x)r^(x) 


So at x = 0 , 


RCN = r (O)r' (0) 
o o 
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Evaluating the above integral for a rounded nose gives 


f 


2 = | r(Or’tf):tf - A + e)~ 2 d£ 


O'' o 


p r o (A)r o 


(A) a - a + e)- 2 


A* 

= RCN (z + e) -2 d£ 

J o 

= rcn(a - -f-} 

ye & + zj 

= RCN 
e 

or e = RCN/2 . 

f 

For a pointed nose, r^(A) = S is finite (S is slope at nose). This gives. 


f 


2 = I r (Or'AOa - A + e)" 2 d? 


o o 


or 


[r;(A)] z J; a ~ A) (? - A + e )" 2 d? 
A 


= S 2 I z(z + e) -2 d£ 

J b 

= S 2 In (5, + e) - lne - 
= -S 2 lne 


-2/S 2 


£ 4- e 


e = e 

So, for a pointed nose, e is very small for a moderate slope. 


Correction Factor K 

For moderately thick bodies, slender body theory underpredicts the maximum thick- 
ness. Figure 2(b) shows the calculated maximum thickness for ellipsoids. For bodies 
that are symmetric fore and aft, that is, for ellipsoids, a factor of K is all that 
is needed to match the maximum thickness. As a result, for ellipsoids, the following 
expression for ^ was used: 
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where 


ip(x,r) 



f Kr (?)r;(?)(x - OdZ 

J AA [(x - C) 2 + r 2 ] l/2 


V, 


AA = A + RCN/ 2 
BB = B - RCT/2 


For an ellipse, RCN = RCT = t 2 /a, where t is the maximum thickness and a is 
the length of the semimajor axis (see fig. 1(b)). In this study, the semimajor axis 
will always be 1/2 and the body will span from 0 to 1. Now, a value of K is 
desired such that r = t at x = 1/2. Setting the above equation for equal to 0 
and expressing the limits of integration in terms of R = 2t 2 , the following equation 
for K is obtained: 


K 


-1 + R + (1 + R 2 ) 1 ^ 2 
. 1 - R + (1 + R 2 ) 1 ^ 2 _ 


+ (1 - R) (1 + R 2 ) 



An approximate form for K(t) can be obtained by expanding the log term and square 
root term. This results in 

f 

K = [1 - R + R ln(R/2) + (3/2) R 2 + 0(R 3 )] -1 


A plot of K versus t is shown in figure 3. Figure 4 shows the amount of over- 
shoot of the stagnation point, using the modified slender body theory. Comparing 
figures 2(a) and 4, the modified theory shows a reduction in the overshoot. 


Compressed Coordinate System 

For the case of more general bodies that are not symmetric fore and aft, a 
compressed coordinate system was introduced to maintain net zero source strength. 
Now, 


q(0 = .tfq SB (£*) 

where 


r* = £ - RCN/2 
^ 1 - RCA 

RCN = radius of curvature of nose 
RCT = radius of curvature of tail 
RCA = (1/2) (RCN + RCT) 


Note: 
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and 


5* = 1 at 
0 < < 1 



q SB 


a*) = u o a * a*) 


Applying this compressed coordinate system to an ellipsoid gives 

q SB^ ^ 2 

q(5) = ~ |^At 2 



1 - RCA q SB (5 ' ) 

= Kq SB (5) 

In the above equation, K is that value calculated previously for an ellipsoid. For 
a more general shape, then, 

q(0 =Mq QB a*) = K(1 - RCA)q SB (5*) 


For these general shapes, K can be calculated based on an ellipsoid with R = RCA. 
Alternatively,^ can be evaluated directly from the requirement that the maximum 
thickness of the body be matched. The following expression for jp was thus used: 


^(x,r) = U t 


<r 2 - 


BB 

*AA 


K(1 - RCA)q(5*)(x - g)dt 
t(x - O 2 + r 2 ] l/2 


Numerical Integration of Equations 

A computer program was written to integrate numerically the equations for the 
stream function and velocities using Simpson’s rule. For example, if a body stream- 
line is desired, the equation for is set equal to zero; then for each value of 
x, a value of r is found by means of an interval-halving scheme such that ip = 0 . 
The body shapes are kept in separate subroutines so that different shapes may be used 
with relative ease. One hundred integration steps (N = 100) was found to be suffi- 
cient in most cases for integrating ip (or the velocities). For locating stagnation 
points, N > 100 was sometimes needed. For body shapes requiring a compressed coor- 
dinate system, £ was transformed to 5* before entering the body-shape subroutines. 
A flowchart of the procedure is shown in figure 5. The subroutine VEL which computes 
the velocity components, is similar to the subroutine STREAM, except that the equa- 
tions for U x and U r are integrated instead of \p. 
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Source Strength Distribution for a Pointed Tail 

Results obtained for the zero streamline using the above expression for \p show 
that the nose is well modeled, but a pointed tail is not. Regular slender body 
theory, however, models the tail shape quite well, but breaks down near the nose. 
Multiplying the source strength by a factor of K(1 - RCA) alone for a body with a 
pointed tail increases the predicted body slope. An additional condition was there- 
fore needed to improve the mathematical model; that is, the derivative of the singu- 
larity strength at a pointed nose or tail was maintained in order to match the body 
slope. A new formulation for q was thus developed for the pointed-tail' case. 


Tail Correction Procedure 
The following equation for the source distribution, 

q(0 = .*q SB (5*) 

models the nose of the body well. But using this distribution gives for the tail 

q* (1) = .*s*’q^ B (l) 

Since slender-body theory models the tail region well, the desired form of q at the 
tail is 

q'd) = q£ B (l) 

Thus, let q take on the following form: 

q(5) = KNq SB a*) + f (?*)(?* - 1) (1 - KN?*')q^ B (l) 


where 

q SB «*) - 2»U o r o (5*)r; (5 *) 

q’ B (l) =2„U o [r;(l)) 2 

and KN is a factor such that r = t at the point of maximum thickness. Note that 
q f ( 1 ) = qsg(l) if f (1 ) = 1/g* 1 ., which is the desired result. For a closed body, 

jf q(5)d5 = 0 

o 

which gives 

J fa*) a* - i)d?* = o 
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In addition, the condition of f (0) = 0 must also be met to maintain the source 
distribution at the nose. 


In summary, f must satisfy the following three requirements: 

1. f(0) = 0 

2. f (1) = 1/?*’ 

3. J fa*)a* - Ddg* = o 

o 


In addition to the above requirements, it was desired that f be a continuous func- 
tion to prevent difficulties in the numerical integration scheme. Therefore, f was 
chosen to be piecewise linear for simplicity. Knowing the desired form of f, the 
following distribution of f was chosen: 


fa*) 



for 0 < K* < ? x 

for ?! < ?* S ? 2 

for ? 2 < ?* < 1 


A sketch of f versus ?* is shown in figure 6. As shown there, there are three 
parameters to be determined: a G , , and £ 2 • By applying condition 3 from above. 


a Q is chosen to be 


3(5, - - 1) - (5* + 5, + -.-.J 

5,(5, + 5,-3) 

For a given body shape, ? 1 and ? 2 are determined by trial and error. Given ? x and 
? 2 , the value of KN can be found by numerically integrating (also using Simpson’s 
rule) the following expression: 

2 j* BB f (?*)(?* - l)(x Q - ?)q’ B (l)d? 

Jaa t(x„ - ?) 2 + t 2 ] l/2 


[q SB (?*)(x o - ?) - f (?*)(?* - l)(x o - ?)q^ B (l)?*']d? 
[ (x„ - ?) 2 + t 2 ] l/2 


Here, x Q is the x-coordinate of the point of maximum thickness. 

The tail correction procedure may be difficult to generalize to other body 
shapes, especially those with abrupt changes in body slopes. Several attempts were 
necessary to choose the and Z 2 values for the cases in this report. 
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Treatment of Test Module Shapes 

The EE and the RTA test modules have blunt ends; however, the radii of 
curvature of the tails are very small relative to the body lengths (see table 1). 

To apply the modified slender body theory with tail correction, the tails of the 
EE and RTA were extended to a point. The extension increased the length of the EE 
by about 7.7%; the RTA length increased by 9.9%. In the analysis thus far, all 
body lengths were normalized to one, as wqs the f ree-stream velocity. Because the 
extensions caused the lengths of the EE and RTA to be greater than one, a "stretch- 
ing" of the coordinate system for these two cases was necessary. A brief explanation 
of this procedure follows. 

Before entering the body shape subroutine, the coordinates are compressed, as 
discussed earlier. The amount by which the EE (or RTA) is extended is computed in 
the body shape subroutine. Before it is entered into the equations that describe the 
body shape, the variable g* is multiplied by the amount 1 of the extension, XE. For 
example, XE for the EE is 1.077308; for the RTA, XE = 1.098950. Also, r 0 (£*XE) 
is multiplied by 1 /XE. The body slope must remain the same as the physical body, so 
rg(£;*XE) is not multiplied by 1/XE. The body streamline generated by this procedure 
has a range of 0 < x < 1 and has a maximum thickness of t/XE, where t is the 
maximum thickness as defined earlier for the actual EE (or RTA) with un extended tail. 
A diagram of the procedure is shown in figure 5(c). Note that 2'iiU 0 r ( i ) r ) I ) and 
27rU 0 [r ) I ) (l) ] 2 are equivalent to q SB and qgg(l) , respectively. ' ’ 


Body Shapes 

Table 1 gives the equations for .the actual body shapes used in this study. 

Table 2 gives the values of , £ 2 , a 0 , arid KN, which produced body streamlines that 
closely matched the actual body shapes. 


Exact Solutions 


Exact solutions describing the flow about a sphere or an ellipsoid can be found 
in the fluid mechanics literature. For a sphere, the potential and stream function 
are given by (see ref. 16) 


4>(x,r) = U x 1 - 


ip(x, r) = U Trr 2 1 + 


£ ’ 2 (x 2 + r 2 ) S/2 ] 
£ + (x 2 + r 2 )’ /2 ] 


where a is the radius of a sphere centered at (x,r) = (0,0). This gives for the 
velocities 


U = U 

X ° (x 2 


U a 3 r 
+ 0 r 2 ) 3/2 L 



U 

y 



*y 

(x 2 + r 2 ) 


5/2 
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u = 

z 


3a 3 xz 

2 / 2 , 2\® A 

(x + r ) 


For a coordinate system centered at the body nose, x should be replaced by (x - a). 
For a sphere of unit length, a = 1/2. 

For flow about an ellipsoid, the potential is given in elliptic coordinates by 
the following expression (see ref. 16): 


<j>(y,v,w) = U Q kyv + Ay 


1 .. (v ± i\ . 

2 V ln yv - lj 1 


where (for v > 1) 


x = f (y , v ,oj) = kjrv 
y = g(y,v,o>) = k cos to (1 - y z )" / ^(v 2 - 1) 
z = h(y,v,oo) **> k sin w(l - y 2 ) 1 ^ 2 (v 2 - 1) 
tan oj = z/y 
k 


2 \ 1 / 2 /-,,2 n 1 / 2 


(a 2 - b 2 ) l/2 


U a 
o 


- h -(Hi) 


1 - e 2 


e = eccentricity 


(a 2 - b 2 ) 


2 x 1/2 


semimajor axis 


A 


b = semiminor axis 
Io obtain the velocity components, the following relations are used: 

u — ^ ^y _i_ 34 1 
x “ 3y 9x 3v 3x 

u — ^ jy . 3<t> 3v 

y 3y 3y 3v 3y 

U = .M. In + jhL 

z 3y 3z 3v 3z 

% 

tote: 3<j)/3dj = 0. From the expression for <J> (y,v) , one obtains 


= U vk + A 
3y o 


|| = U Q yk + Ay 


& u ln (N - l ] 

[i -H - H-r] 
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In order to find the partial derivatives with respect to the Cartesian coordi- 
nates, the following matrix inversion is necessary: 
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Inverting the matrix results in 








v (1- 

y 2 ) 

9y a 
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. 

-y sin 03(1 - y 2 ) 1/ ^ 2 (v 2 • 
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- y 2 ) 
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l) 
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V COS 03 
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v sin 03(1 - y 2 ) 1 / 2 ( v 2 - 

9x 

k(v 2 - 
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k(v 2 - y 2 ) 

9(jo _ 
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9 03 _ 



-sin 03 





9 03 

COS 0) 

9x 
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9y " 

k(i 
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- 1) 1/2 


9z 

k(l - y 2 ) l/z (v 2 - 1) 1/2 


Hence, the velocity components are given by 

U = U + v Ft in / -- ±-j \ - ~~ 1 

x o k I 2 yv - 1 J v 2 - 1 1 

A /I 2 \ ^ 

y _ -Ay cos 03(1 - y ) 
y k(v 2 - y 2 )(v 2 - l) 1 ^ 2 

U = "Ay sin 03(1 - y 2 ) 1 ^ 2 
Z k(v 2 - y 2 ) (v 2 - l) l/2 

Given (x,y,z) , v can be found, which then gives y, as follows 

y 2 + z 2 = k 2 (1 - y 2 )(v 2 - 1) 

= k 2 (v 2 + y 2 - y 2 v 2 - 1) 

= k 2 v 2 + x 2 /v 2 - x 2 - k 2 
or 

o - »* - v t g+y2 k ^ 2+k2 ) + (!)* 

The above quadratic equation results in four roots. Since v must be greater than 1, 
three of the roots are eliminated and so: 
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- j x 2 + y 2 + z 2 + k 2 
2k 2 


x 2 + y 2 + z 2 + k 2 
^ k 2 


) -u] 1 


/ 2 


The value of k is determined by the maximum thickness of the ellipse (b = t); the 
value of the semimajor axis will always be 1/2. 


Incorporation of Body-Induced Velocities into a 
Comprehensive Rotor Analysis 

After the various body models were developed to a satisfactory degree, the next 
step was to introduce the body-induced velocities into the comprehensive helicopter 
analysis CAMRAD (ref. 17). The solution procedure for the rotor behavior is basi- 
cally unchanged, since no iteration between the body and rotor is involved (except to 
update the body-induced velocities when the rotor position relative to the body 
changes as the trim iteration proceeds) , The cases considered in this report involve 
a single rotor, and a body at zero angle of attack. Hence, the wind-tunnel configu- 
ration of reference 17 can be used, with zero pitch and yaw angles of the body axes. 
The axis systems used are shown in figure 7. The calculation of the body-induced 
velocities Is performed for each location on the rotor disk. The blade position -en 
the rotor disk (scaled with the body length) relative to the body nose is (see 
fig. 7): 

r = ^ -r B + r R + R SF r tP I 


where 

£/R = body length/rotor radius 

r^ - position of the body origin relative to the body c.g. (body axes) 

D 7 

r^ = position of the rotor hub relative to the body c.g. (body axes) 

r^ - blade position relative to the hub location (shaft axes) 

Rg-p = matrix which transforms body-axis system (origin at c.g.) to shaft-axis system 
(origin at hub position); defined by shaft pitch angle 0 r and cant angle <j) r 

The calculation of is defined in reference 17 (p. 143). The blade position r is 

sent to a subroutine that calculates the body-induced velocities. The velocities are 
then- included as perturbations in the rotor section velocity components. Since 
CAMRAD computes the rotor section velocities relative to the shaft axes, it is neces- 
sary to transform the body-induced velocities to the shaft axes system. Hence, 

* = r sf^ (or) 


where 

V = (V x ,Vy,V z ) = body-induced velocity relative to body axes, 
scaled by free-stream velocity 
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SIR - rotor tip speed 

A - induced velocity at rotor disk relative to shaft axes 

The location on the rotor disk at which the body-induced velocities are calcu- 
lated depends on the blade motion; therefore, ideally the induced velocities should 
be calculated at least once per circulation iteration within CAMRAD. However, this 
resulted in a large number of induced velocity calculations which were found to be 
unnecessary in achieving convergence. Calculation of the induced velocities once per 
control iteration was found to be more than sufficient, since the rotor position 
changes little during the final stages of the trim iteration, as shown in figure 8, 
for a representative forward flight analysis. Options to compute the velocities at 
other locations in CAMRAD are possible. 

Figures 9(a) and 9(b) show the entry points of the body-induced velocity sub- 
routine into CAMRAD. To assess the influence of the various body shapes discussed 
thus far, eight options were programmed. These eight options and the necessary 
inputs required by CAMRAD are shown in figure 9(c). A list of the subroutines in 
CAMRAD, which were modified, and a list of new subroutines, which were added, are 
shown in figure 9(d), A skeleton of the subroutine BODY I is shown in figure 9(e). 
Complete details, of the CAMRAD program structure and input are given in reference 18. 

For the case of OPSHAP = -1 or -2, an unformatted file of velocities is read; 
no calculations of body-induced velocities are made. The velocity file is generated 
by a panel code or by some other means independent of CAMRAD. The disadvantage of 
this method for accounting for body-induced effects is that the velocities cannot be 
updated if the rotor blade position relative to the body changes during the analysis. 
The file has the logical name VELOCITYFILE, and will be read by a statement of the 
form 


READ (NFILE) ( ( (VINT(K,I , J) ,K = 1,3), I = 1,MRA),J = 1 ,MPSI) 

where K is the index for the three velocity components, I is the index for the 
blade radial station, and J is the index for the blade azimuth position. The 
velocity should be calculated at the midpoints of the blade aerodynamic panels defined 
by the parameter RAE in namelist NLRTR (see ref. 18); the total number of stations, 
MRA, is also from namelist NLRTR. The azimuth positions are ip = 360 * J/MPSI (deg), 
where the_total number MPSI is from namelist NLTRIM. For OPSHAP = -1, a velocity 
file of V(r,\J;) relative to the (Xp,Yp,Zp) coordinate system is read (see fig. 7). 

The velocities are scaled by the free-stream velocity. These velocities will be 
transformed to A as described above. For OPSHAP = -2 3 a velocity file of A(r,^) 
relative to the (Xg,Y s ,Zg) coordinate system is read. In this case, the velocities 
are scaled by the tip speed QR. 


ANALYTICAL MODEL OF ROTOR 


Three types of rotors (table 3) were used in the present study. Nonuniform 
inflow and a free wake geometry (for two or three revolutions of the wake) were used 
to model the rotor environment. Results obtained from use of a prescribed wake 
geometry showed no significant change from the free wake results. Aerodynamic calcu- 
lations were performed at 15 radial stations along the blade and at every 15° of 
azimuth. The radial stations along the blade become more concentrated near the blade 
tip. For rotor A, a teetering flap mode with one harmonic of motion was used. For 
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rotors B and C, six flap/lag bending modes (including the rigid body modes) and one 
(blade) pitch mode were used; seven harmonics of motion were used for each mode. 

Rotor A is a stiff model rotor, for which only the performance was calculated. Both 
rotor performance and loads were calculated for rotors B and C. 

All the calculations were performed for sea level standard conditions. The 
rotor tip Mach number was held constant (see table 3), and the advance ratio was in 
the range p =0.15 to 0.50, depending on the case considered (see table 4). The 
rotor Jhrust and the tip-path-plane tilt relative to the rotor shaft were trimmed to 
specified values by varying the rotor collective pitch and cyclic pitch control 
angles. The thrust was varied from C^/o = 0. 05 to 0.09, typically, and the lateral 
flapping angle was always trimmed to zero. The longitudinal flapping angle (and 
hence the tip-path-plane angle) were trimmed to values depending on the case consid- 
ered (see table 4) . 

The body velocities were recalculated after each control change during the rotor 
trim iteration, ensuring that the velocities were evaluated at the correct, that is, 
at the final trimmed rotor position. 

The adequacy of the rotor model can be assessed by comparing the theory with 
experimental rotor data. Figure 10 shows the L/D as a function of thrust for rotor A 
over a range of advance ratios. The data were obtained from reference 6. The corre- 
lation is good, although the experimental data were not accurate enough to show a 
trend with tip-path-plane angle (for a fixed advance ratio). Figure 11 shows L/D 
versus advance ratio for two values of C^/a for rotor B. Figure 12 shows the one- 
half peak-to-peak blade loads for rotor B. Figures 12(b) and 12(c) were calculated at 
the 70% and 60% blade radial station, respectively. The data were obtained from 
reference 19. Both the data and the theory for figures 11 and 12 include the effects 
of the RTA, but as will be shown later in the paper, the RTA has a very small effect 
on the rotor behavior. The correlation is good to fair; similar correlation is fehown 
in reference 11. 

No data are available yet for the isolated behavior of rotor C. 


Nonuniform Inflow 

The changes in blade loads and profile power are due to changes in the complex 
flow field of the rotor blade. The perturbation velocities are fixed by the body. 
These velocities do not simply produce an incremental change in the blade loads and 
power. Because of the nonlinear relationship between the perturbation velocities and 
the blade loads and the profile power, proper calculation of the blade angle of 
attack distribution is necessary. Changes in the induced power caused by the body 
would be zero if a uniform inflow model (momentum theory with empirical corrections) 
was used. Even the direct interference power requires a detailed thrust distribution 
over the rotor disk. Hence, a nonuniform inflow model is necessary to assess body- 
induced effects on rotor behavior. 


Wake Geometry 

As stated earlier, a free wake geometry was used for all calculations presented 
in this paper. Comparisons were made with prescribed wake results, and no significant 
changes were found in performance or loads caused by the wake distortion. It follows 
that body-induced wake geometry variations would not be important either. Therefore, 
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the absence of the distortion of the wake geometry by the body in the present analy- 
sis is not significant. The change in the body flow field due to the rotor, which 
then changes the perturbation velocities at the rotor, is also neglected. 


Trimming (Convergence of Solution) 

Results presented here are for a fixed tip-path-plane angle, a^pp . An alternate 
approach would be to trim to a specified propulsive force. The presence of the body 
does change the propulsive force for a fixed a T pp; therefore, slightly different 
results are obtained for trimming to propulsive force rather than trimming to a«ppp. 
These differences were not found to be significant for the cases presented here, 
which are for small propulsive force. 


RESULTS: VERIFICATION OF APPROACH 


The accuracy of the modified slender body theory will be examined for case 1 of 
table 4: ellipsoids and axisymmetric bodies with a NACA four-digit airfoil thickness 

distribution. The ellipsoid case is useful because the exact solution for the body- 
induced velocities is known. The airfoil shapes are a reasonable approximation for 
typical wind tunnel test modules. The complete modified slender body theory with 
tail correction will serve as the "exact 11 solution for the airfoil shaped bodies. 
Since the purpose is to establish the limits of the method, high values of the body 
thickness (up to t = 0.50) are used. 


Ellipsoids 


The validity of the models developed using slender body theory can be determined 
by first considering the flow field of ellipsoids, for which the exact solution is 
available. 

Figure 13 shows the dividing streamline (body shape) of ellipsoids, as obtained 
by unmodified and modified slender body theory. Compared with the exact body shape, 
the modified theory shows good accuracy for the 60%-thick body, and fair accuracy for 
the 80%-thick body. The 100%-thick body is even less accurate, but the modified 
theory always matches the maximum thickness. 

Figures 14 and 15 show the flow field (velocity magnitude and angle of attack, 
respectively) about a 60%-thick ellipsoid in a plane above the body. Figures 16 
and 17 and figures 18 and 19 show the flow field about 80%- and 1 100%-thick ellipsoids, 
respectively. The body-plane separation for all three bodies is 10% of the body 
length. Relative to the exact solution, slender body theory shows significant errors 
for all three bodies. The modified theory predicts both the velocity magnitude and 
angle of attack well, even though the body shape was not matched exactly (fig. 13). 
Note that for the 60%-thick body, the modified theory gives results that are nearly 
identical to those of the exact solution. 

Figures 20-25 show the section angle of attack change induced by the body on a 
rotor blade for 60%-, 80%-,, and 100%-thick ellipsoids. The rotor radius is equal to 
the body length, and the rotor-body separation is equal to 10% of the body length. 

The advance ratio is 0.4. The rotor hub is located above the body center in fig- 
ures 20, 22, and 24; it is located above the body quarter chord in figures 21, 23, 
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and 25, Again, the modified slender body theory results are almost identical to 
those of the exact solution for 60%- and 80%-thick bodies. The modified theory 
results for the 100%-thick body show fair agreement with the exact solution. The 
body-induced changes in the velocity magnitude seen by the rotor are not significant, 
as shown by figure 26 (note that fig. 26(d) is for the isolated rotor). Here, the 
body is 80%-thick, the advance ratio is 0.4, and the rotor hub is located above the 
body center. The parameters for figures 27 and 28 are the same as in figures 22 
and 23 except that the advance ratio is 0.2. At the lower advance ratio, the modi- 
fied slender body theory results are also quite accurate. 


Airfoil Shaped Bodies 

The applicability of the theoretical method will be examined by comparing the 
predicted and actual dividing streamlines for more general (and more useful) body 
shapes than the ellipsoids. Typical results are as follows. Figure 29 shows axi- 
symmetric bodies with a NACA four-digit airfoil thickness distribution. The 10%-thick 
body is well modeled by slender body theory except for the nose stagnation point 
location, which is correctly placed by using the modified theory. For the 30%- and 
5.0%-thick bodies, results are shown from the modified theory with and without the 
tail correction.' It is seen that the full, modified slender body theory gives good 
results. , * 

The angle of attack change at the rotor blade is shown in figure 30 for a 
30%-thick airfoil. The advance ratio is 0.4, the body length is equal to the rotor 
radius, and the rotor-body separation is again equal to 10% of the body length. The 
rotor is located above the body maximum thickness point for the symmetrical airfoil; 
this point is located at the body 0. 3-chord station, Aft of the hub position, the 
full modified theory is actually closer to slender body theory than it is to the 
modified theory without the tail correction (because slender body theory models the 
shape of the tail correctly). Forward of the hub position, the modified theory 
without the tail correction and slender body theory solutions are quite similar. 


Error in Performance and Loads 

After satisfactory modeling of the ellipsoids and airfoils had been achieved, 
the next step was to determine the difference in rotor performance and loads as cal- 
culated using modified slender body theory and the exact solution (for airfoil shapes, 
modified slender body theory with tail correction results served as the "exact" 
solution). All values in this paper for edgewise and flapwise bending moments are 
one-half peak-to-peak values taken from steady-state time histories. 


Figures 31(a) and 31(b) show, as a function of thrust, the error in the calcu- 
lated lateral cyclic control angle and the rotor lift-to-drag ratio, respectively, 
caused by ellipsoids of various thicknesses. Rotor A was used for these calculations; 
other parameters are given in table 4, case 1. The errors are given by 


A(e lc ) 


■ < v - 


( 6 ) 

ic exact body 


4(L/D) - (L/D) - (L/D) exact body 


IS 



Clearly, slender body theory shows significant errors, except for thin ellipsoids. 
The modified slender body theory results are good for moderately thick bodies. 
Figure 32 shows similar results for airfoil shaped bodies. The errors are given by 

A ^ 9 ic^ ^ic^ ^lc^MSBT with tail correction 


A(L/D) = (L/D) - (L/D) 


MSBT with tail correction 


Figures 33(a) and 33(b) are plots of the errors in the oscillatory edgewise and 
flapwise bending moments, respectively, as a function of thrust for various ellip- 
soids. The moments were calculated at the 60% radial station. Rotor B was used for 
these calculations together with the values given in case 1 of table 4. The errors 
were calculated using the same procedure for calculating the performance parameters 
that was used above and are plotted as a fraction of the oscillatory loads for the 
exact body solution. Even for large thickness ratios, modified slender body theory 
shows very little error compared with the exact solution for both edgewise and flap- 
wise bending moments. Slender body theory, however, shows large errors for both 
cases. Figure 34 shows similar results for the edgewise bending moment for various 
airfoil shapes. The flapwise bending moment showed little effect of the body, so no 
error calculations are presented. For a 50%-thick airfoil, modified slender body 
theory shows an improvement over slender body theory. Without the tail correction, 
there is an error of about 7.5% for the edgewise bending loads. For the 10%- and 
30%-thick airfoils, slender body theory gives errors that are less than 5%. 


Influence of Body on Rotor 

The effects of various bodies on rotor performance and loads were next analyzed. 
Ideal body shapes are considered here since the results provide insight into the 
mechanisms involved in the body-induced changes of rotor performance and loads. 

Figures 35(a) and 35(b) show the changes in lateral cyclic control and lift-to-drag 
ratio, respectively , caused by various body shapes. Rotor A was used together with 
the parameters in case 1 of table 4. The ellipsoid results were calculated using the 
exact solution; the airfoil shaped body results were calculated using the modified 
slender body theory with tail correction. Figure 35(a) shows that for a given body 
shape, there is relatively no change in A8 lc with rCtor thrust. 

The body produces a basic once-per-revolution 'variation of the vertical inter- 
ference velocity seen by the rotor blade, with maximum amplitude at azimuth angles 
ip = 0° and 180°. Hence a lateral cyclic pitch change is needed to maintain the 
tip-path-plane lateral inclination at zero. Figure 35(b) shows the decrease in L/D 
caused by the body. Figure 36 shows the lift-to-drag ratio of rotor A caused by 
various body shapes. The performance loss increases in magnitude as the body thick- 
ness increases. The loss is significant only for large thickness ratios. 

Rotor B, together with the parameters of case 1, was used to generate the results 
shown in figures 37 and 38. Figure 37 shows the changes in lift-to-drag ratio, 
oscillatory edgewise bending load, and oscillatory flapwise bending load for various 
ellipsoids and airfoil shapes. Unlike rotor A, rotor B shows an increase in perfor- 
mance up to about C«j7a - 0.075 (compare figs. 37(a) and 35(b)). In figure 37, the 
blade oscillatory bending moments are plotted as a fraction of the isolated rotor 
loads. The bending loads increase for all body shapes over the given range of thrust. 
Figure 38 shows the edgewise and flapwise oscillatory loads as a function of thrust. 
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Pitch link loads were also calculated, but no significant influence of the body 
was found (5%— 10% change for the ellipsoids, 3%-5% change for the airfoil shapes). 
Also, there was little difference in the pitch link loads predicted using the various 
body theories compared above. 


Source of Performance Loss Caused by Body 

To determine the cause of the performance change produced by the body, a detailed 
analysis of the effect of an 80%- thick ellipsoid on rotors A and B was performed. 

Table 5 shows the breakdown of Cp/o (power coefficient divided by solidity) for 
rotor A as affected by an 80%-thick ellipsoid. The rotor thrust was Cp/o = 0.08 
(case 1). From table 5, it follows that there is a small increase in the performance 
directly caused by the interference velocity. The loss caused by the profile power 
and induced power is larger, however, than the gain caused by the interference power. 
Thus, the lift-to-drag ratio decreases. Note that the same behavior is shown for 
both tip-path-plane angle trim and propulsive force trim. v 

Polar plots were generated for X-£ n t (interference inflow ratio = V^ nt /fiR), 
a (blade section angle of attack), c^ (drag coefficient), Cp Q (profile power coeffi- 
cient), c^ (lift coefficient), Cf (thrust coefficient), C Pi (induce^ power coeffi- 
cient), and Cp^ nt (interference power coefficient) to determine how each parameter 
varied over the rotor disk. A brief summary of the plots is discussed below. 

The vertical interference velocity X.£ nt was found to be negative (up) on the 
front of the disk, and positive (down) on the rear of the disk. Hence, Xi n t was 
roughly proportional to cos \p (a 1/rev variation with for/aft symmetry). Lateral 
cyclic, A0 1q , is needed to cancel the resulting change in flap moment, AMp, resulting 
from this 1/rev variation. The lateral cyclic control angle, A0 1C , and Xint are 
both proportional to cos rp, but 0 is independent of r and X is larger near the 
hub. Roughly, X^ nt = (A t /r 2 )cos tJj, where At is the value of A at the tip. The 
flap moment is given by 

AM F ~ r 2 (r0 - X^) = (r 3 A0 lc - X t )cos ip 

It follows that the total blade flap moment change, which must be zero for an articu- 
lated rotor, is 



Now, the change in angle of attack can be expressed as 

Aa = 0 — — = X fc cos \p ^ 
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Thus , we expect that Aa = 0 at r = 4 = 0,63. Results show Aa to be zero at 

approximately r = 0.75. 

In general, a was found to increase inboard and decrease outboard on the front 
of the disk. On the rear of the rotor disk, a decreased inboard and increased out- 
board. Figure 39 shows the angle of attack change caused by the body over the rotor 
disk for this case. This angle of attack change produces a change in the blade sec- 
tion drag coefficient, but the pattern is different because c<j is a nonlinear 
function of a. The drag coefficient c^ was found to increase inboard of the region 
of ip = 180°. On the advancing side, however, c^ increased outboard. The change 
in a caused by the body on the retreating side of the disk shifts the stall region 
by a small amount. The lower values of a on the rear of the disk reduce the drag 
somewhat. Although a is negative in the tip region in the second quadrant of the 
rotor disk, the magnitude is such that the cj is increased on most of the advancing 
side in the tip region. The profile power coefficient Cp Q was found to vary simi- 
larly to ,Cj. 

The rotor-blade lift coefficient c% and the thrust coefficient C«p were found 
to vary similarly to a. The maximum circulation (which determines the wake-^induced 
velocities) is found to be greater at both ip = 0° and at ip - 180° at constant 
C T /a , so the radial distribution of the circulation is being changed. The induced 
velocity was somewhat higher on the front and back of the disk when the body was 
included. f 


The induced power coefficient Cp_^ was found to vary similarly to C«p. There 
is a small net increase in Cp^ since generally increases from the front to the 

back of the rotor disk. The interference power coefficient Cp^ varies similarly 


to Ai nt . The rotor thrust coefficient C-p is higher on the front of the disk 
(partly because of body effects), so the power increase from the front of the disk 
is greater than the power loss from the rear of the disk. 


Table 6 shows the breakdown of C p /a for rotor B as affected by an 80%-thick 
ellipsoid (case 1). The same basic behavior in rotor A is also observed in rotor B. 
A high value of C j/o is needed, however, in order for the sum of the increases in 
induced and profile power to be greater than the interference power (note that 
C y/o = 0.08 for table 5). 


Source of Bending Load Increase Caused by the Body 

To analyze the body effects on the rotor loads, rotor B was used in combination 
with an 80%-thick 1 ellipsoid . Figures 40(a) and 40(b) show the time histories of the 
edgewise and flapwise bending loads, respectively. Figure 40 was generated under 
the conditions of case 1, with C -p/a = 0.08 (see also table 6) . The bending moments 
were calculated at the 60% radial station. Figure 40 shows that the effect of the 
body is mainly on the front of the disk from ip = 90° to ^ = 270°. Referring to 
figures 37(b) and 37(c), the increase in the edgewise and flapwise bending moments 
caused by the body is 74% and 26%, respectively, for an 80% ellipsoid at Cj/o = 0.08. 
The change in the edgewise bending moment is positive on the front of the disk, corre- 
sponding to a lead motion. The change in the flapwise bending moment is negative in 
this region, corresponding to a downward flapping motion. As stated earlier, a was 
found to increase inboard and decrease outboard on the front of the disk, as did C T 
and c^, which explains the behavior of the bending moments. 


21 



Figures 41(a) and 41(b) show the modulus of the harmonics of the oscillatory 
edgewise and flapwise bending moments, respectively, as a function of the number of 
the harmonic, n. The 1/rev is the dominant harmonic. There is an increase in the 
modulus of both bending moments because of the body for n = 1 to n =3 . The peak 
in the area of n = 5 for the edgewise bending moment is due to a coupled flap-lag 
blade natural frequency near 5fl. The effect of the body is thus not just a local 
phenomenon, but rather changes the entire pattern of blade loading. 


RESULTS: ROTOR-BODY AERODYNAMIC INTERFERENCE 


The effects of the body on rotor performance and loads will be examined for 
cases 2 to 5 of table 4. Even though the modified slender body theory has been 
developed here only (for axisymmetrie bodies at zero angle of attack, a number of 
practical configurations can be analyzed: a typical full-scale rotor test in the 

Aines 40- by 80-Foot Wind Tunnel; the Ames small-scale interactional aerodynamics 
tests; the future Ames full-scale interactional aerodynamics tests; and a rotor- 
hybrid airihip configuration. 


Wind Tunnel Test Modules 

i ’ 

The next step in the analysis was to model the test modules used in the 40- by 
80-Foot Wind Tunnel at Ames and to determine their effects on rotor performance* and 
loads . 

Figure 42 shows the shapes of the RTA and EE test modules. Modified slender 
body theory modeled these shapes with no significant error in body streamlines (recall 
that these test modules are modeled with the tails extended to a point) . All the 
following results for the influence of the RTA and EE test modules were obtained 
using the modified slender body theory with tail correction. 

The surface pressure distribution on a 1 /6-scale EE, obtained from reference 8, 
is compared with the distribution calculated by the modified slender body theory with 
tail correction in figure 43. The theory compares well with the data and with the 
panel code from reference 8, except for x > 0.95. An explanation for the lack of 
agreement in this vicinity is that the data represent the occurrence of flow separa- 
tion. The theory, at this stage of development, does not model separation. The 
influence of the rotor on the body pressure is greater than the error caused by 
separation at the tail (see ref. 8). 


Full-Scale Test Configuration # 

The combination of rotor B and the RTA is a configuration typical of a full- 
scale wind tunnel test in the Ames 40- by 80-Foot Wind Tunnel. Such tests are 
intended to measure the isolated rotor characteristics. Figures 44 to 50 were gener- 
ated under the conditions of case 2 in table 4. Figure 44 shows the change in angle 
of attack produced by the RTA for the baseline values of case 2. The advance ratio 
is 0.4. 

Figure 45 shows the increase in the lift-to-drag ratio of rotor B caused by the 
RTA for various advance ratios. As shown, there is a negligible performance improve- 
ment caused by the RTA\(A (L/D) = 0 to 0.2), The profile and induced power are basically 
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unaffected by the presence of the body; the change in L/D is due almost entirely to 
the interference power. 

Figures 46(a) and 46(b) show the increase in the oscillatory edgewise and flap- 
wise bending moments, respectively. The edgewise moments, which were calculated at 
the 60% radial station, show a 10%-15% increase because of the body. The flapwise 
bending moments, which were calculated at the 70% radial station, show a 5%-10% 
increase caused by the RTA. The calculated influence of the RTA on the pitch- link 
loads was negligible (±2%). 

Figure 47 shows the lift-to-drag ratio of rotor B with and without the presence 
of the RTA for various advance ratios. As noted earlier, very little change is caused 
by the body. Figures 48(a) and 48(b) show thh edgewise and flapwise bending moments, 
respectively. 

Figures 49(a) and 49(b) show the time histories of the edgewise and flapwise 
bending moments, respectively. The advance ratio is 0.5 and the Cj/a is 0.07. The 
edgewise moment shows small changes at a few azimuth locations. The flapwise moment 
shows almost no change as the blade progresses around the rotor disk. Figure 50 
shows the moduli of the harmonics of the edgewise bending moment as a function of 
the harmonic number. 


Small-Scale Interactional Aerodynamic Test Configurations 

To simulate some of the small-scale rotor-body interactional tests conducted at 
Ames (refs. 6-8), the effects of scaled-models of the RTA and EE on rotor A were 
analyzed. 

Figures 51(a) and 51(b) show the change in angle of attack of the rotor blade 
caused by the EE for the two longitudinal hub positions in case 3. The baseline 
values for the advance ratio and the rotor-body vertical separation were used. 

Moving the hub position aft, as shown in figure 51(b), changes the angle of attack 
distribution significantly. 

The baseline conditions of case 3 in table 4 apply to table 7 and to figures 51 

to 55. The parameter variations for case 3 were chosen to match those in reference 7. 

The baseline configurations for the RTA and EE approximately model the full-scale 
test configurations. 

Figure 52 shows the lift-to-drag ratio for rotor A as affected by the RTA and EE 
for the baseline configuration at two different advance ratios. The decrease in L/D 
caused by the EE increases as the advance ratio is increased. The RTA causes a 
small increase in L/D at the higher advance ratio. 

Figure 53 shows the effect of changing hub position on the lift-to-drag ratio. 
The body used is the EE. Moving the hub aft with respect to the body nose creates a 
positive increase in L/D as compared with the isolated rotor. The advance ratio 
for figure 53 is 0.3. Table 7 shows the contributions to the power at C T /a =0.07 

for the isolated rotor case, and for the cases of rotor A in combination with the RTA 

and EE. There are small increases in the profile and induced power, but the change 
in the total performance is due primarily to the direct interference power. 

Figures 54(a) and 54(b) show the change in L/D caused by variations in advance 
ratio, rotor height, and hub position for the RTA and EE, respectively. In 
figure 54(a), the variation in rotor height has little effect on the change in L/D. 
Moving the rotor aft with respect to the body nose increases the change in L/D. The 
increase in advance ratio from p = 0.15 to 0.2 causes a decrease in A(L/D). As 
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y is increased from 0,2 to 0.3, however, A (L/D) is increased. For the EE, increas- 
ing the rotor height causes a small increase in A (L/D) as shown in figure 54(b). 

As the hub position is moved aft, there is an increase of A (L/D) of about 0.4. The 
increase in advance ratio causes A (L/D) to decrease. 

Figure 55 shows the lift-to-drag ratio for rotor A as affected by the EE for 
various tip-path-plane angles at an advance ratio of 0.3. As a T pp becomes more 
negative, L/D is seen to decrease. Figures 56(a) and 56(b) show the change in L/D 
as affected by the RTA and EE, respectively, for several tip-path-plane angles at an 
advance ratio of 0.3. Figure 56(a) does not show any discernible trend of A (L/D) 
with a T pp, whereas in figure 56(b), A(L/D) is seen to become less negative as c^pp 
becomes more negative. 

Although these calculations are for the small-scale test configurations of 
references 6-8, it is not possible to validate the theory using the rotor performance 
data of reference 6, a result of the systematic error now believed to exist in the 
data. Reference 6 gives data for the rotor alone and with the RTA model. Figure 4 
of reference 6 shows no trend of the isolated rotor performance with a-ppp, whereas 
the calculations show A (L/D) =0.5 to 1.0 for an 8° change in c^PP (fig. 55). 
Figures 5(a) and 8(b) of reference 6 show changes in A (L/D) of the order of 0.5 for 
the range of advance ratios, hub positions, and rotor-body separation distances inves- 
tigated; but the calculations show an influence of these parameters thpt is much less 
(fig. 54(a)). Hence, all that can be concluded is that the calculations support 
the conclusion that there are systematic errors of the order of 0.5 in the measured 
rotor L/D. 


Full-Scale Interactional Aerodynamic Test Configurations 

The combination of EE and rotor C is a configuration to be used as part of the 
full-scale interactional aerodynamic tests at Ames (case 4 of table 4). Table 8 
shows the contributions to the power coefficient for the baseline values of case 4. 
Comparing tables 7 and 8, shifting the hub position aft has a less favorable effect 
on the interference power for rotor C than for rotor A. Also, the induced power 
shows an increase as a result of the aft shift in hub position for rotor C, whereas 
rotor A shows a decrease in the induced power as well as in the profile power. 

Table 8 shows an overall net performance gain of 2% caused by the aft shift in hub 
position, whereas table 7 shows a 6% net performance gain. 

Figure 57 shows the effect of the hub position on the lift-to-drag ratio at two 
advance ratios. Figure 58 shows the effect of hub position on the change in L/D 
caused by the presence of the EE at two advance ratios. The effect of the EE is 
seen to be very small (A (L/D) = ±0.1). Figures 59(a) and 59(b) show the increase 
caused by the EE in the oscillatory edgewise and flapwise bending moments, respec- 
tively, as affected by the hub position at two advance ratios. The moments were 
calculated at the 50% radial station. The change in the edgewise bending moment 
decreases as the advance ratio is increased (hub position fixed) , and increases as 
the hub position is moved aft (advance ratio fixed). The change in the flapwise 
bending moment, however, increases as the advance ratio is increased. 

Figures 60(a) and 60(b) show the oscillatory edgewise and flapwise bending 
moments, respectively, as affected by the hub position at two advance ratios. The 
edgewise moment shows an increase of 10%-20% caused by the body for the baseline hub 
position and a 25%-45% increase for the aft hub position. The flapwise moment shows 
an increase of 5%-15% because of the body. There is little effect of the hub position 
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on the flapwise moment. The calculated effect of the EE on the pitch link loads 
was small (a 4%-8% increase). 

Figures 61 and 62 show the time histories of the bending moments at the 50% 
radial station for C -p/cr = 0.07. The advance ratio is 0.3 and 0.4 for figures 61 
and 62, respectively. The effect of the EE on the edgewise moment is greatest at 
the front of the rotor disk. The flapwise moment shows small changes caused by the 
body at various locations on the rotor disk. Figure 63 shows the moduli of the 
harmonics of the edgewise moment as a function of the harmonic number. 


Rotor-Airship Interference 

Because an airship is similar in shape to an ellipsoid, performance calculations 
were made for a typical hybrid airship configuration using a 5 0%- thick ellipsoid. 

Case 5 of table 4 lists the parameter values used for the calculations. This air- 
craft operates at low speed (p = 0.15 here). Four rotors, two forward and two aft, 
were placed at the bottom. edge of the body. Figure 64 shows the lift-to-drag ratio 
of a front rotor and a rear rotor. The presence of such a large body produces a sig- 
nificant effect on the L/D. The L/D of the rear rotor is seen to increase by 
approximately 17%, and the L/D of the front rotor decreases by 15%. The effect of 
the body had negligible influence on the profile and the induced power. The effect 
on power required was due almost entirely to the interference power.' Note that the 
induced velocity is downward at the front of the body and upward at the rear of the 
body. 


COMPUTATION EFFICIENCY 


The computational efficiency of modified slender body theory can be assessed by 
comparing the time required for the body analysis and the time needed for the rotor 
portion of the analysis. The amount of time required for the body analyses depended 
on the frequency at which the body velocities were updated within CAMRAD. Two 
approaches were used to estimate the computational efficiency of the body analysis. 

1. Update body velocities (over entire rotor disk) once per wake iteration 
(0PINTV =2). 

2. Update body velocities once per control change in trim iteration (0PINTV = 3 

One wake iteration involves typically 5 to 15 steps of control changes (see fig. 8), 
Note that the rotor analysis begins with a uniform inflow solution, then proceeds to 
a nonuniform inflow with prescribed wake geometry (iterated if desired), and then to 
a nonuniform inflow with free wake geometry (iterated if desired) (see fig. 9(a) 
and refs. 17-18). 

Approach (2) can be considered a fully converged body velocity calculation, 
since the rotor motion, hence the rotor position relative to the body, changes very 
little during the final steps of the trim iteration (see fig. 8). The following 
table shows the percentage of the time required for the entire job that was used in 
the wake influence coefficient calculations and the body induced velocity calcula- 
tions using (1) and (2). Approach (2) requires more updates (a factor of 10-20) of 
the body-induced velocities than approach (1). 
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Approach 

Wake influence 
coefficient 
calculation, % 

Body velocity 
calculation, % 

i 

35-40 

5-10 

2 

15-25 

40-70 


Regarding the accuracy of approaches (1) and (2), approach (1) with two wake 
iterations gave less than a 1% difference in performance and loads than approach (2). 
Also, approach (1), with two iterations, required about the same total time, or even 
less, than approach (2), with one wake iteration. For an optimally accurate and 
computationally efficient solution, a more complex method of determining when to 
update the body-induced velocities is needed. It would be sufficient to update the 
body velocities three to five times at the beginning of the trim iteration. Updating 
the velocities at the end of each trim iteration degrades the efficiency signifi- 
cantly, as seen by the results of approach (2). Listed below is a summary of rela- 
tive computation times required by the different analyses* 


Analysis 

Relative time 

Body-induced velocities calculated 
by modified slender body theory 

15-25 

Free wake geometry 

60-100 

Wake influence coefficient matrix 

150-200 

Complete rotor solution (each wake 
iteration) 

400-600 

Body-induced velocities calculated 
by panel method 

1000-3000 


If the body velocities are updated 3-5 times per wake iteration in the rotor 
solution, then the body velocities, as calculated by modified slender body theory, 
will require 10%-20% of the time required for a complete rotor solution. If 
a panel method is used, however, the time required to calculate the velocities will 
be 2-6 times that' required for a complete rotor solution, depending on the complexity 
of the panel code. (Although a more efficient panel code could be constructed for 
the case of an axisymmetric body at zero angle of attack, the above comparison of 
computation time refers to the more general problem.) 


CONCLUSIONS 


Modified slender body theory was used to assess body-induced effects on rotor 
behavior. The effects of several body shapes on three different rotors were studied. 
The wide range of configurations covered was made possible by the computational effi- 
ciency of the body model. Conclusions drawn from this study are listed below. 

1. Moderately thick ellipsoids are modeled well by modified slender body theory. 
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2. The streamlines of body shapes with pointed ends can be closely matched by 
the modified slender body theory with tail correction, but a trial and error proce- 
dure is necessary to find and £ 2 , 

3. In assessing the effects of simple body shapes on the rotor, slender body 
theory produced significant errors (relative to the exact solution) except for thin 
shapes, whereas the modified slender body theory was accurate for moderately thick 
bodies* 


4. Since the body- induced effects depend on the detailed aerodynamic environ- 
ment of the rotor blade, nonuniform inflow induced velocity calculations must be used 
in the rotor analysis. A free wake model produced essentially the same results as a 
prescribed wake geometry model, implying that the neglect of body-induced distortions 
of the wake geometry is acceptable. 

5. The calculated influence of the body on the rotor performance was generally 
small for the cases considered here. The direct interference power can be positive 

or negative, depending on the body configuration. The interference power is typically 
not zero, even for cases with exact fore-aft symmetry of the interference velocity, a 
result of the asymmetry of the rotor loading distribution. For the cases considered 
here, the profile power and induced power were always increased by the presence of 
the body. The profile power increase was generally greater than the induced power 
increase. The net performance change caused by the body could be favorable or 
unfavorable, depending on the sign of the direct interference power and its magnitude 
relative to the profile power change. 

6. The calculated oscillatory blade bending moments were always increased by 
the body for the cases considered here. The edgewise bending moment changes were 
larger than the flap wise bending moment changes. In some cases, the loads were 
increased by a significant fraction of the isolated rotor loads. The interference 
effects were due to the general changes in the rotor flow pattern, rather than to 
localized effects of the body. The influence of the body on the blade loads occurred 
primarily on the front of the rotor disk. No significant effects on pitch link 
loads were found. 

7. For the case of a typical full-scale rotor test in the 40- by 80-Foot Wind 
Tunnel, negligible effects- of the Rotor Test Apparatus (RTA) on the rotor performance 
and pitch link loads were calculated. Only small increases in the blade bending 
moments were found. Hence, such a test does produce essentially isolated rotor 
behavior. 

8. For the case of the Ames small-scale aerodynamic interference tests, a per- 
formance change of about A(L/D) - 0.2 was predicted for the Easter Egg (EE) body 
shape. A 6% performance change (A(L/D) = 0.4) was calculated for the two body longi- 
tudinal positions tested. The corresponding influence of the RTA body shape was 
small. The calculated influence of rotor-body vertical separation was small for both 
body shapes. 

9. For the case of the Ames full-scale aerodynamic interference tests, it is 
predicted that the EE body will produce negligible performance changes. The oscilla- 
tory edgewise bending moments will be 10%-20 % higher than for the isolated rotor with 
the EE at the baseline position. If the EE body were shifted forward, the loads 
would be 25%-45% higher than for the isolated rotor. Smaller effects are predicted 
for the influence of the body on the flapwise bending moment. 
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10* For the case of a rotor-hybrid airship, a significant performance increase 
and decrease were predicted for the rear and front rotors, respectively (A (L/D) = ±15%). 

11. The computational efficiency of the body model is one or two orders of 
magnitude better than that of a typical potential flow panel code. The body analysis 
using the modified slender body theory required only 10% to 20% of the computation 
time required for the rotor analysis. 


RECOMMENDATIONS 


The predicted influence of the body on the rotor should be verified by compari- 
son with experimental data. Such comparisons should be done initially for simple 
shapes, such as the RTA and EE wind tunnel test modules, for which the accuracy of 
the isolated body analysis is well established. 

The modified slender body theory should be extended to eliminate the restrictions 
to axi symmetric bodies and zero body angle of attack. It is anticipated that the 
computational efficiency of the method will be easily maintained, but that some loss 
in accuracy representing the body streamline will result. The present results imply 
that for practical shapes the accuracy will remain acceptable. Techniques for 
modeling flow separation should be established, and a more general approach for the 
pointed-tail correction would be desirable. The incorporation of rotor-induced 
effects on the body aerodynamics should also be considered for the modified slender 
body theory. 
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TABLE 1.- EQUATIONS FOR VARIOUS BODY SHAPES 
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TABLE 2.- VALUES OF Ci , Z z > a 0 , AND KN FOR DIFFERENT BODY SHAPES 


Body 

Thickness 
ratio, % 

t 

RCN 

RCT 



a o 

KN 

Extended 

length 

Ellipsoid 

60 

0.30 

0. 18 

0. 18 

. 


— 

2.3227 


Ellipsoid 

80 

.40 

.32 

.32 

— 

— 

— 

4.5377 

— 

Sphere 

100 

.50 

.50 

.50 

— 


— 

12.8526 

— 

NACA airfoil 

10 

.05 

.0110 

0 

0.10 

0.95 

0.4629 

1.2680 

— 

NACA airfoil 

30 

.15 

.0992 

0 

. 10 

.95 

.4629 

2.6539 

— 

NACA airfoil 

50 

.25 

.2755 

0 

. 10 

.95 

.4629 

5.5858 

— 

RTA 


.08745 

.0585 

0 

.05 

. 06 

15.4504 

.9261 

1.0990 

EE 


. 1537 

.6774 

0 

. 10 

.65 

.8530 

2.2642 

1.0773 


TABLE 3.- ROTOR TYPES USED IN ANALYSIS 


Type 

Rotor A 
(teetering) 

Rotor B 
(articulated) 

Rotor C 
(articulated) 

Number of blades 

2 

4 

4 

Radius, m 

1.12 

6.71 

7.01 

Solidity ratio 

.0651 

. 0748 

.0705 

Twist, deg 

-10 

-10 

(nonlinear) 

-14 

(nonlinear) 

Lock number 

3.44 

9.08 

7.19 

Tip Mach number 

.60 

.60 

.70 
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TABLE 5.- PERFORMANCE LOSSES FOR ROTOR A AT C /a = 0.08 WITH 80% -THICK ELLIPSOID (CASE 1) 
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TABLE 6.- PERFORMANCE CHANGES FOR ROTOR B WITH 80%-THICK ELLIPSOID (CASE I) 
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TABLE 7.- PERFORMANCE CHANGES FOR ROTOR A WITH 
RTA AND EE (CASE 3) 


Trim to C^/o = 0.07 



c Pi /a 

(induced power) 
C Pint /a 

(interference power) 

c p 0 /a 

(profile power) 


(C D , + C t 


Cp/a = 


J Pi _ Pint 

a 




A(C p ./a) 

4 (Cpmt/a) 

A (Cp^/a ) 
ACp / o 
L/D 


Isolated 

rotor 

With 

RTA 

With 

EE 

With 

EE a 

0.00102 

0.00103 

0.00106 

0.00104 

— 

-.00005 

... 00001 

-.00012 

.00178 

.00181 : 

.00182 

. 00181 

.00280 

.00279 

.00289 

.00273 

— 

.00001 

.00004 

.00002 

— 

-.00005 

.00001 

-.00012 

— 

.00003 

. 00004 

.00003 

— 

-.00001 

.00009 

-.00007 

7.50 

7.54 

7.26 

7.70 


0.471 
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TABLE 8.- PERFORMANCE CHANGES FOR ROTOR C 
WITH EE (CASE 4) 



Trim to C-j/a = 

0.07 


Isolated 

rotor 

With 

EE 

With 

EE a 

c Pi /a 

0. 00116 

0.00120 

0.00122 

(induced power) 




C Pint /a 

— 

-.00004 

-.00011 

(interference power) 




Cp 0 /a 

.00185 

.00187 

,00187 

(profile power) 

(C p± + Cp int + Cp o ) 

s /a = 

.00301 

,00304 

.00298 

A(C p ./a) 

— 

.00004 

. 00006 

A(C Plnt /a) 

— 

-.00004 

-.00011 

A (0 p /o) 

: 

.00002 

. 00002 

A (C p /a) 

— 

.00003 

-.00003 

L/D 

6.99 

6.92 

7.05 


a (x/£) hub = 0.500. 
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LOOKING AFT 


(b) 

Figure 1.- Coordinate system and geometry used to define body shapes, 
(a) General axisymmetric body; (b) prolate spheroid. 
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(a) Amount of overshoot of stagnation point as a percent of body length using 

slender body theory for ellipsoids. 



(b) Calculated maximum thickness using slender body theory for ellipsoids. 
Figure 2.- Amount of overshoot of stagnation point and calculated maximum thickness. 


39 






Figure 4 



- Amount of overshoot of stagnation point as a percent of body length 
using modified slender body theory for ellipsoids. 



J 
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MAIN PROGRAM 


INITIALIZE 

N 

U o 

t 

x 

r 1 

i p 

Y o 

CALL STREAM (x, r v \p^) 

5 r 2 = r-j - Ar 

CALL STREAM (x r r 2 , ^ 2 ) 

IF (V^-! < ^ 0 < ^2> G0 T0 10 
r 1 =r 2 

V'l = ^2 

GO TO 5 

10 rr = 0 

15 r = (r-j + r 2 )/2 

CALL STREAM (x, r, 4>) 

IF ( |rr - r| < t) GOTO 100 

IF (\p <^ Q ) r-j = r 
IF {\p >*P Q ) r 2 = r 
rr ~ r 

GO TO 15 

100 PRINT x.r.tf/ 

CALL VEL (x, r, U x , U r ) 

PRINT U x , U r 

STOP 

END 

(a) Summary of main program. 

Figure 5 . — Procedure for calculating streamlines and velocities. 
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SUBROUTINE STREAM (x f r, p) 


1 1 _ BB AA 
H 

a = o 

DO 1 00 I = 1 , N 
>.. g-RCN/2 

5 1 - RCA , - 

SUBROUTINE BODY SHAPE (£*, q) 

_____ 

r o = r o 

q = Vo* 

RETURN 

100 % = P (x, r, £, q ({*)) 

P = Z \pft 

RETURN 

(b) Summary of stream function and body shape subroutines. 

Figure 5.- Continued. 
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MAIN PROGRAM 

SUBROUTINE STREAM (x, r. \p) 

• 

E ~ RCn 

? 2 

1 - RCA 

SUBROUTINE BODY SHAPE (£*) 

I COMPUTE XE 

COMPUTE f (l-j, ?2' KN already known at this stage) 

t* = xE,r 

r o = r o<£*)/XE 
r o = r o (£*> 

r = ?*/xe 

q = KN*r 0 + f,(£* - 1).(1 - KN,$*')*(r 0 '{1)) 2 
1 RETURN 

CALCULATE $ 

RETURN 


END 


(c) Procedure for inclusion of an extended tail. 
Figure 5.- Concluded. 
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Figure 6.- Function f for tail correction 


(X, Y, 2): COORDINATE SYSTEM WITH ORIGIN AT BODY NOSE 

(X s , Y $ , Z s ): COORDINATE SYSTEM WITH ORIGIN AT HUB 

(X F , Yp, Zp): COORDINATE SYSTEM WITH ORIGIN AT CG 



Figure 7,- Blade position on rotor disk relative to body origin 




UNIFORM INFLOW | NONUNIFORM INFLOW 





PRESCRIBED WAKE 

FREE 

GEOMETRY 

■ WAKE 



0 5 10 15 20 25 30 


CONTROL ITERATION NUMBER 

Figure 8.- Convergence of rotor blade position during typical trim iteration. 
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TRIM 


TIMER 

IF TRIM RESTART, GO TO RESTART ENTRY POINT 
I F OPINTV = 1, CALL BODY! 

UNIFORM INFLOW 
IF ITERU^O 

IF OPINTV = 2, CALL BODY! 

TRIMS 

NONUNIFORM INFLOW, PRESCRIBED WAKE GEOMETRY 
FOR IT = 1 TO ITERR 

WAKEC1 

WAKEC2 

IF OPINTV = 2, CALL BODY! 

TRIMI 

NONUNIFORM INFLOW, FREE WAKE GEOMETRY 
FOR IT = 1 TO ITERF 

WAKEC1 

WAKEC2 

IF OPINTV = 2, CALL BODY! 

TRIMI 

• 


(a) BODY I subroutine calls in TRIM., 
Figure 9 . - CAMRAD for program modifications. 
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RAMF 


TIMER 

IF OPINTV = 3, CALL BODY! 

FOR COUNTC “ 1 TO ITERC (CIRCULATION ITERATION) 
IF OPINTV = 4, CALL BODY I 


FOR COUNTM - 1 TO ITERM (MOTION ITERATION) 
IF OPINTV = 5, CALL BODY I 

* 

TEST MOTION CONVERGENCE 
TEST CIRCULATION CONVERGENCE 

BODYF 

TIMER 

■(b) BODYI subroutine calls in RAMF. 
Figure 9.- Continued. 
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PARAMETER 

DESCRIPTION 

FSBODY 

BLBODY 

WLBODY 

(x, y, z) COORDINATES OF BODY NOSE RELATIVE TO THE HUB POSITION, 
NORMALIZED BY R 

LENGTH 

BODY LENGTH, 2/R 

THICK 

BODY THICKNESS RATIO, 2t/8 (FOR OPSHAP = 1-6 ONLY) 

OPSHAP 

INTEGER SPECIFYING BODY SHAPE: 

1 - ELLIPSOID (EXACT SOLUTION) 

2 - ELLIPSOID (MSBT) 

3- ELLIPSOID (SBT) 

4 - NACA OOxx (MSBT WITH TAIL CORRECTION) 

5 - NACA OOxx (MSBT) 

6 - NACA OOxx (SBT) 

7 - ROTOR TEST APPARATUS (MSBT WITH TAIL CORRECTION) 

8 - EASTER EGG (MSBT WITH TAIL CORRECTION) 

-1 - READ VELOCITY FILE (RELATIVE TO BODY AXES, SCALED WITH V) 

-2 - READ VELOCITY FILE (RELATIVE TO SHAFT AXES, SCALED WITH J2R) 

NMSBT 

NUMBER OF NUMERICAL INTEGRATION STEPS USED TO INTEGRATE 
EQUATION FOR BODY-INDUCED VELOCITY 

OPINTV 

INTEGER SPECIFYING WHEN BODY-INDUCED VELOCITIES 
ARE UPDATED: 

0 - NOT CALCULATED AT ALL (NO BODV) 

1 - ONCE PER CASE 

2 - ONCE PER WAKE ITERATION 

3 - ONCE PER CONTROL ITERATION 

4 - ONCE PER CIRCULATION ITERATION 

5 - ONCE PER MOTION ITERATION 


(c) New input parameters in namelist NLBODY and new common block INTAER (blade 
position calculation also uses OPWKBP parameter from namelist NLWAKE). 

Figure 9.- Continued. 


49 


















MODIFIED SUBROUTINES 


SUBROUTINE 

MODIFICATIONS 

INPTB 

INCLUSION OF NEW PARAMETERS IN NAMELIST NLBODY 

PRNTB 

NEW INPUT PARAMETERS PRINTED 

TRIM 

CALLS BODYI 

RAMF 

CALLS BODYI 

LOADR1 

POLAR PLOT INCREMENTS 

TIMER 

NEW TIME SUMMARIES 


NEW SUBROUTINES 


SUBROUTINE 

DESCRIPTION 

BODYI 

CALCULATES BODY-INDUCED VELOCITIES OVER ROTOR DISK 

BODYB1 

CALCULATES BLADE POSITION (IN SHAFT AXES) AT AERODYNAMIC 
STATIONS; DERIVED FROM SUBROUTINE WAKEB1 

VBODY 

CALCULATES BODY-INDUCED VELOCITY 

EXACT 

CALCULATES FLOW ABOUT ELLIPSOID - EXACT SOLUTION 

VMSBT 

INTEGRATES VELOCITY EQUATION 

CURVE 

CALCULATES BODY NOSE/TAIL RADIUS OF CURVATURE 

OPTION 

SELECTS BODY SHAPE 

ELL 

ELLIPSOID- USES MSBT 

ELLSB 

ELLIPSOID- USED SBT 

NACAMSBT 

NACA OOxx— USES MSBT WITH TAIL CORRECTION 

KVALUE 

CALCULATES VALUE OF KN FOR NACA OOxx 

FZETA 

CALCULATES PARTS OF INTEGRAND FOR INTEGRATION IN KVALUE 

NACAMSB 

NACA OOxx— USES MSBT 

NACASB 

NACA OOxx— USES SBT 

RTA 

ROTOR TEST APPARATUS - USES MSBT WITH TAIL CORRECTION 

EGG 

EASTER EGG - USES MSBT WITH TAIL CORRECTION 


(d) List of new and modified subroutines. 
Figure 9. - Continued. 
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BODY I 


IF OPSHAF > 0 
BODYB1 


VBODY 


VMSBt 

CURVE 


OPTION 


ELL 

ELLSB 


NACAMSBT 

KVALUE 

FZETA 


NACAMSB 


NACASB 

RTA 

EGG 


EXACT 


I F OPSH AP < 0 

READ VELOCITY FILE 


(e) Skeleton of BODYI subroutine (called by TRIM and RAMF) . 

Figure 9.- Concluded. 
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FORWARD 

SPEED 



Figure 14.- Plot of velocity magnitude divided by free stream velocity in plane above 
body — 60% ellipsoid. (a) Exact theory; (b) slender body theory; (c) modified 
slender body theory. 
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Figure 16.- Plot of velocity magnitude divided by free stream velocity in plane above 
body — 80%-thick ellipsoid. (a) Exact theory; (b) slender body theory; (c) modi- 
fied slender body theory. 
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Figure 20.- Plot of rotor blade angle of attack change (deg) due to a 60%-thick 
ellipsoid in plane of rotor disk — hub position above body mid-chord, p = 0,4-. 
(a) Exact theory; (b) slender body theory; (c) modified slender body theory. 
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Figure 21.- Plot of rotor blade angle of attack change (deg) due to a 60%-thick 
ellipsoid in plane of rotor disk — hub position above body quarter-chord, \x : 
(a) Exact theory; (b) slender body theory; (c) modified slender body theory. 












Figure 22.- Plot of rotor blade angle of attack change (deg) due to an 80%-thick 
ellipsoid in plane of rotor disk — hub position above body mid-chord, y = 0,4, 
(a) Exact theory; (b) slender body theory; (c) modified slender body theory. 
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Figure 23.- Plot of rotor blade angle of attack change (deg) due to an 80%-thick 
ellipsoid in plane of rotor disk — hub position above body quarter-chord, y = 
(a) Exact theory; (b) slender body theory; (c) modified slender body theory. 
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Figure 24,- Plot of rotor blade angle of attack change (deg) due to a 100%-thick 
ellipsoid in plane of rotor disk — hub position above body mid-chord, p =0.4. 
(a) Exact theory; (b) slender body theory; (c) modified slender body theory. 
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Figure 25.- Plot of rotor blade angle of attack change (deg) due to a 100%-thick 
ellipsoid in plane of rotor disk — hub position above body quarter-chord, y = 
(a) Exact theory; (b) slender body theory; (c) modified slender body theory. 































Y 

- 1.0 -.5 0 .5 1.0 



V 

- 1.0 -.5 0 .5 1.0 



<b) 



figure 27.- Plot of rotor blade angle of attack change (deg) due to an 80%-thick 
ellipsoid in plane of rotor disk — hub position above body mid-chord, y = 0.2. 
(a) Exact theory; (b) slender body theory; (c) modified slender body theory. 
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RADIUS RADIUS RADIUS 


EXACT BODY SHAPE 

SLENDER BODY THEORY 

MODIFIED SLENDER BODY THEORY 

MODIFIED SLENDER BODY THEORY 

(WITHOUT TAIL CORRECTION) 



-.2 0 .2 .4 .6 .8 t.O 1.2 

(c) NACA 0050 CHORD 


Figure 29.- Body streamlines for axisymmetric bodies with NACA 4-digit airfoil 
thickness distribution. (a) NACA 0010; (b) NACA 0030; (c) NACA 0050. 
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Figure 30.- Plot of rotor blade angle of attack change (deg) due to a 30%-thick 
airfoil shaped body in plane of rotor disk, y = 0.4. (a) Slender body theory; 

(b) modified slender body theory; (c) modified slender body theory with tail 
correction. 
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(b) Lift-to-drag ratio. 
Fi gure 31 . - Cone luded . 
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SLENDER BODY THEORY 

modified slender body theory 

(WITHOUT TAIL CORRECTION) 


THICKNESS 

RATIO 



C-y/o 

(a) Lateral cyclic control angle (deg). 

Figure 32.- Error in lateral cyclic control angle (deg) and lift-to-drag ratio for 
rotor A in combination with airfoil shaped bodies of various thicknesses (case 1) . 


74 



SLENDER BODY THEORY 
MODIFIED SLENDER BODY THEORY 
(WITHOUT TAIL CORRECTION) 




(b) Lift-to-drag ratio. 
Figure 32.- Concluded . • 
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SLENDER BODY THEORY 

MODIFIED SLENDER BODY THEORY 



(b) Oscillatory f lapwise bending loads. 


Figure 33.- Concluded. 
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(b) Lift-to-drag ratio. 
Figure 35.- Concluded. 
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G-j-/(7 


(c) Oscillatory f lapwise bending loads (at 0.6R), 
Figure 37.- Concluded. 
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.0010 r 


.0008 H 


.0006 h 


ISOLATED ROTOR 

WITH ELLIPSOID 

WITH AIRFOIL SHAPE 

THICKNESS 

RATIO 



.0004 


.0002 h 


n I I I : I I I 

.04 .05 .06 .07 .08 ' .09 

Cj/a 


(b) Oscillatory flapwise bending loads 
Figure 38.- Concluded. 
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Figure 42.- Body streamlines of test modules (dashed lines are tail extensions) 
(a) Rotor Test Apparatus (RTA) ; (b) Easter Egg (EE) . 














80 

deg 


360 


moment at 0.7R. 

ts for rotor B and RTA configuration at 
= 0.07 (case 2). 
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Figure 55.- Lift-to-drag ratio of rotor A ai 

path-plane angles (case 3). 
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Figure 59 Effect of hub position on the increase in oscillatory bending moments at 
0.5R for rotor C due to the EE for two advance ratios (case 4). (a) Oscillatory 

edgewise bending moment; (b) oscillatory flapwise bending moment. 







.0008 


.0007 


.0006 


.0005 


E * .0004 
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.0002 


.0001 



ISOLATED ROTOR 

( x /C)|_|ub = °- 500 

— — — - (x/e) HUB = 0.206 


/ 


/ 



H = 0.3 



Cy /<; 

(b) Oscillatory flapwj.se bending moment . 
Figure 60.- Concluded. 
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Figure 64.- Lif t-to-drag ratio of a front and rear rotor of a typical airship 
modeled as a 5 0%- thick ellipsoid (case 5). 
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